Multiscale Frozen Density Embedding/Molecular Mechanics Approach for Simulating Magnetic Response Properties of Solvated Systems.

We present a three-layer hybrid quantum mechanical/quantum embedding/molecular mechanics approach for calculating nuclear magnetic resonance (NMR) shieldings and J-couplings of molecular systems in solution. The model is based on the frozen density embedding (FDE) and polarizable fluctuating charges (FQ) and fluctuating dipoles (FQFµ) force fields and permits the accurate ab initio description of short-range nonelectrostatic interactions by means of the FDE shell and cost-effective treatment of long-range electrostatic interactions through the polarizable force field FQ(Fµ). Our approach's accuracy and potential are demonstrated by studying NMR spectra of Brooker's merocyanine in aqueous and nonaqueous solutions.

UV-Resonance Raman Spectra of Systems in Complex Environments: A Multiscale Modeling Applied to Doxorubicin Intercalated into DNA.

UV-Resonance Raman (RR) spectroscopy is a valuable tool to study the binding of drugs to biomolecular receptors. The extraction of information at the molecular level from experimental RR spectra is made much easier and more complete thanks to the use of computational approaches, specifically tuned to deal with the complexity of the supramolecular system. In this paper, we propose a protocol to simulate RR spectra of complex systems at different levels of sophistication, by exploiting a quantum mechanics/molecular mechanics (QM/MM) approach. The approach is challenged to investigate RR spectra of a widely used chemotherapy drug, doxorubicin (DOX) intercalated into a DNA double strand. The computed results show good agreement with experimental data, thus confirming the reliability of the computational protocol.

Amide Spectral Fingerprints are Hydrogen Bonding-Mediated.

The origin of the peculiar amide spectral features of proteins in aqueous solution is investigated, by exploiting a combined theoretical and experimental approach to study UV Resonance Raman (RR) spectra of peptide molecular models, namely N-acetylglycine-N-methylamide (NAGMA) and N-acetylalanine-N-methylamide (NALMA). UVRR spectra are recorded by tuning Synchrotron Radiation at several excitation wavelengths and modeled by using a recently developed multiscale protocol based on a polarizable QM/MM approach. Thanks to the unparalleled agreement between theory and experiment, we demonstrate that specific hydrogen bond interactions, which dominate hydration dynamics around these solutes, play a crucial role in the selective enhancement of amide signals. These results further argue the capability of vibrational spectroscopy methods as valuable tools for refined structural analysis of peptides and proteins in aqueous solution.

An efficient and robust procedure to calculate absorption spectra of aqueous charged species applied to NO2.

Accurate calculation of absorption spectra of aqueous NO2- requires rigorously sampling the quantum potential energy surfaces for microsolvation of NO2- with at least five explicit water molecules and embedding the resulting clusters in a continuum solvent accounting for the statistical weighted contributions of individual isomers. This method, which we address as ASCEC + PCM, introduces several desired features when compared against MD simulations derived QM/MM spectra: comparatively fewer explicit solvent molecules to be treated with expensive QM methods, the identification of equilibrium structures in the quantum PES to be used in further vibrational spectroscopy, and the unequivocal identification of cluster orbitals undergoing electronic transitions and charge transfer that originate the spectral bands.

A polarizable three-layer frozen density embedding/molecular mechanics approach.

We present a novel multilayer polarizable embedding approach in which the system is divided into three portions, two of which are treated using density functional theory and their interaction is based on frozen density embedding (FDE) theory, and both also mutually interact with a polarizable classical layer described using an atomistic model based on fluctuating charges (FQ). The efficacy of the model is demonstrated by extending the formalism to linear response properties and applying it to the simulation of the excitation energies of organic molecules in aqueous solution, where the solute and the first solvation shell are treated using FDE, while the rest of the solvent is modeled using FQ charges.

Multilevel Density Functional Theory.

Following recent developments in multilevel embedding methods, we introduce a novel density matrix-based multilevel approach within the framework of density functional theory (DFT). In this multilevel DFT, the system is partitioned in an active and an inactive fragment, and all interactions are retained between the two parts. The decomposition of the total system is performed upon the density matrix. The orthogonality between the two parts is maintained by solving the Kohn-Sham equations in the MO basis for the active part only, while keeping the inactive density matrix frozen. This results in the reduction of computational cost. We outline the theory and implementation and discuss the differences and similarities with state-of-the-art DFT embedding methods. We present applications to aqueous solutions of methyloxirane and glycidol.

Binding of SARS-CoV-2 to Cell Receptors: A Tale of Molecular Evolution.

The magnified infectious power of the SARS-CoV-2 virus compared to its precursor SARS-CoV is intimately linked to an enhanced ability in the mutated virus to find available hydrogen-bond sites in the host cells. This characteristic is acquired during virus evolution because of the selective pressure exerted at the molecular level. We pinpoint the specific residue (in the virus) to residue (in the cell) contacts during the initial recognition and binding and show that the virus⋅⋅⋅cell interaction is mainly due to an extensive network of hydrogen bonds and to a large surface of noncovalent interactions. In addition to the formal quantum characterization of bonding interactions, computation of absorption spectra for the specific virus⋅⋅⋅cell interacting residues yields significant shifts of Δλmax =47 and 66 nm in the wavelength for maximum absorption in the complex with respect to the isolated host and virus, respectively.

Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems.

Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.

Calculation of Linear and Non-linear Electric Response Properties of Systems in Aqueous Solution: A Polarizable Quantum/Classical Approach with Quantum Repulsion Effects.

We present a computational study of polarizabilities and hyperpolarizabilities of organic molecules in aqueous solutions, focusing on solute-water interactions and the way they affect a molecule's linear and non-linear electric response properties. We employ a polarizable quantum mechanics/molecular mechanics (QM/MM) computational model that treats the solute at the QM level while the solvent is treated classically using a force field that includes polarizable charges and dipoles, which dynamically respond to the solute's quantum-mechanical electron density. Quantum confinement effects are also treated by means of a recently implemented method that endows solvent molecules with a parametric electron density, which exerts Pauli repulsion forces upon the solute. By applying the method to a set of aromatic molecules in solution we show that, for both polarizabilities and first hyperpolarizabilities, observed solution values are the result of a delicate balance between electrostatics, hydrogen-bonding, and non-electrostatic solute solvent interactions.

Molecular spectroscopy of aqueous solutions: a theoretical perspective.

Computational spectroscopy is an invaluable tool to both accurately reproduce the spectra of molecular systems and provide a rationalization for the underlying physics. However, the inherent difficulty to accurately model systems in aqueous solutions, owing to water's high polarity and ability to form hydrogen bonds, has severely hampered the development of the field. In this tutorial review we present a technique developed and tested in recent years based on a fully atomistic and polarizable classical modeling of water coupled with a quantum mechanical description of the solute. Thanks to its unparalleled accuracy and versatility, this method can change the perspective of computational and experimental chemists alike.

Vibrational circular dichroism under the quantum magnifying glass: from the electronic flow to the spectroscopic observable.

We present a comprehensive methodology for the analysis and interpretation of vibrational circular dichroism spectra supported by novel graphical representations. The tools rely on the vibrational transition current density (VTCD) associated with a molecular vibration, whose visualization allows exploration of the physical origin of the electronic contribution to the electric and magnetic vibrational dipole transition moments. Different ways of visualizing VTCD from 2D maps to 3D representations are reported and applied to molecular systems of growing complexity. An extension of the VTCD analysis to fully anharmonic spectra within the second-order vibrational perturbation theory (VPT2) is discussed. The analysis is applied to different types of chiral systems: the doubly deuterated oxirane (2S,3S)-oxirane-d2, taken as a reference to validate our implementation; 1,3-difluoroallene, a simple rigid system to explore the application of VTCD to anharmonic VCD spectra. Finally, the analysis of VTCD has been used to better understand the origin of the signal enhancement in peptides linked to the ferrocene groups. The merits and shortcomings of the methods are discussed, and some perspectives for future developments are offered.

A combined experimental and theoretical study of optical rotatory dispersion for (R)-glycidyl methyl ether in aqueous solution.

The dispersive optical activity for aqueous solutions of non-rigid (R)-glycidyl methyl ether (R-GME) has been explored synergistically from experimental and theoretical perspectives. Density functional theory analyses performed with the polarizable continuum model for implicit solvation identified nine low-lying stable conformers that are interconverted by rotation about two large-amplitude torsional coordinates. The antagonistic chiroptical signatures predicted for these structural isomers were averaged under a Boltzmann-weighting ansatz to estimate the behavior expected for a thermally equilibrated ensemble. This led to optical rotatory dispersion profiles that reproduced the overall shape of observations but failed to achieve uniform agreement with measured specific-rotation values even when anharmonic vibrational corrections were applied. A mixed QM/FQ paradigm, whereby quantum-mechanical (QM) calculations of optical activity were combined with classical molecular dynamics simulations of explicit solvation that included mutual-polarization effects by means of fluctuating charges (FQ), was enlisted to elucidate the microsolvation environment and gauge its impact upon conformer distributions and response properties. Although quantitative accord with experiments remained elusive, this approach revealed strong variations in the magnitude and sign of rotatory powers for R-GME as the configuration of surrounding water molecules evolved, thereby highlighting the inherently dynamical nature of the solvated chiroptical response, calling into question the validity of "static" descriptions based on the presumption of distinct energy minima, and giving insight into the inherent complexity posed by the modeling of such properties for solvated systems.

Computational simulation of vibrationally resolved spectra for spin-forbidden transitions.

In this computational study, we illustrate a method for computing phosphorescence and circularly polarized phosphorescence spectra of molecular systems, which takes into account vibronic effects including both Franck-Condon and Herzberg-Teller contributions. The singlet and triplet states involved in the phosphorescent emission are described within the harmonic approximation, and the method fully takes mode-mixing effects into account when evaluating Franck-Condon integrals. Spin-orbit couplings, which are responsible for these otherwise forbidden phenomena, are accounted for by means of a relativistic two-component time-dependent density functional theory method. The model is applied to two types of chiral systems: camphorquinone, a rigid organic system that allows for an extensive benchmark, and some members of a class of iridium complexes. The merits and shortcomings of the methods are discussed, and some perspectives for future developments are offered.

Polarizable Embedding Approach for the Analytical Calculation of Raman and Raman Optical Activity Spectra of Solvated Systems.

We present an analytical formulation and implementation of Raman and Raman Optical Activity (ROA) spectra within a three-layer fully polarizable Quantum Mechanical (QM)/ Molecular Mechanics (MM)/Polarizable Continuum Model (PCM) approach. Polarization effects in the MM layer are modeled by exploiting the Fluctuating Charges (FQ) method, in which MM solvent atoms are endowed with electric charges that can be mutually polarized by the solute QM density. Because of its fully polarizable atomistic description, QM/FQ/PCM is able to account for specific solvent effects like those due to hydrogen bonds, providing a physical picture for protic solvents such as water. Applications to aqueous (R)-methyloxirane and (S)-methyllactate are presented, and results are compared with available experimental data.

Two-Component Noncollinear Time-Dependent Spin Density Functional Theory for Excited State Calculations.

We present a linear response formalism for the description of the electronic excitations of a noncollinear reference defined via Kohn-Sham spin density functional methods. A set of auxiliary variables, defined using the density and noncollinear magnetization density vector, allows the generalization of spin density functional kernels commonly used in collinear DFT to noncollinear cases, including local density, GGA, meta-GGA and hybrid functionals. Working equations and derivations of functional second derivatives with respect to the noncollinear density, required in the linear response noncollinear TDDFT formalism, are presented in this work. This formalism takes all components of the spin magnetization into account independent of the type of reference state (open or closed shell). As a result, the method introduced here is able to afford a nonzero local xc torque on the spin magnetization while still satisfying the zero-torque theorem globally. The formalism is applied to a few test cases using the variational exact-two-component reference including spin-orbit coupling to illustrate the capabilities of the method.

Effective Inclusion of Mechanical and Electrical Anharmonicity in Excited Electronic States: VPT2-TDDFT Route.

We present a reliable and cost-effective procedure for the inclusion of anharmonic effects in excited-state energies and spectroscopic intensities by means of second-order vibrational perturbation theory. This development is made possible thanks to a recent efficient implementation of excited-state analytic Hessians and properties within the time-dependent density functional theory framework. As illustrated in this work, by taking advantage of such algorithmic developments, it is possible to perform calculations of excited-state infrared spectra of medium-large isolated molecular systems, with anharmonicity effects included in both the energy and property surfaces. We also explore the use of this procedure for the inclusion of anharmonic effects in the simulation of vibronic bandshapes of electronic spectra and compare the results with previous, more approximate models.

Relativistic Two-Component Particle-Particle Tamm-Dancoff Approximation.

With the recent introduction of the particle-particle random-phase and Tamm-Dancoff approximations to ab initio theory, routine queries of traditionally difficult systems, such as diradicals and doubly excited states, have been made possible. However, although a wealth of inquiry has been directed to investigating these methods, the current formulations have been restricted to spin-collinear systems, leaving the methods incapable of treating noncollinearity and spin-orbit relativistic effects in excited states. In this work, we extend the particle-particle Tamm-Dancoff approximation to suit two-component Hamiltonians to explicitly treat relativistic effects in excited states. After reviewing the theory and computational implementation, we demonstrate the accuracy of this extension by evaluating the fine structure splittings some of atomic and molecular systems.

Real time propagation of the exact two component time-dependent density functional theory.

We report the development of a real time propagation method for solving the time-dependent relativistic exact two-component density functional theory equations (RT-X2C-TDDFT). The method is fundamentally non-perturbative and may be employed to study nonlinear responses for heavy elements which require a relativistic Hamiltonian. We apply the method to several group 12 atoms as well as heavy-element hydrides, comparing with the extensive theoretical and experimental studies on this system, which demonstrates the correctness of our approach. Because the exact two-component Hamiltonian contains spin-orbit operators, the method is able to describe the non-zero transition moment of otherwise spin-forbidden processes in non-relativistic theory. Furthermore, the two-component approach is more cost effective than the full four-component approach, with similar accuracy. The RT-X2C-TDDFT will be useful in future studies of systems containing heavy elements interacting with strong external fields.

Direct Atomic-Orbital-Based Relativistic Two-Component Linear Response Method for Calculating Excited-State Fine Structures.

In this work, we present a linear-response formalism of the complex two-component Hartree-Fock Hamiltonian that includes relativistic effects within the Douglas-Kroll-Hess and the Exact-Two-Component frameworks. The method includes both scalar and spin relativistic effects in the variational description of electronic ground and excited states, although it neglects the picture-change and explicit spin-orbit contributions arising from the two-electron interaction. An efficient direct formalism of solving the complex two-component response function is also presented in this work. The presence of spin-orbit couplings in the Hamiltonian and the two-component nature of the wave function and Fock operator allows the computation of excited-state zero-field splittings of systems for which relativistic effects are dominated by the one-electron term. Calculated results are compared to experimental reference values to assess the quality of the underlying approximations. The results show that the relativistic two-component linear response methods are able to capture the excited-state zero-field splittings with good agreement with experiments for the systems considered here, with all approximations exhibiting a similar performance. However, the error increases for heavy elements and for states of high orbital angular momentum, suggesting the importance of the two-electron relativistic effect in such situations.

The consequences of improperly describing oscillator strengths beyond the electric dipole approximation.

The interaction between a quantum mechanical system and plane wave light is usually modeled within the electric dipole approximation. This assumes that the intensity of the incident field is constant over the length of the system and transition probabilities are described in terms of the electric dipole transition moment. For short wavelength spectroscopies, such as X-ray absorption, the electric dipole approximation often breaks down. Higher order multipoles are then included to describe transition probabilities. The square of the magnetic dipole and electric quadrupole are often included, but this results in an origin-dependent expression for the oscillator strength. The oscillator strength can be made origin-independent if all terms through the same order in the wave vector are retained. We will show the consequences and potential pitfalls of using either of these two expressions. It is shown that the origin-dependent expression may violate the Thomas-Reiche-Kuhn sum rule and the origin-independent expression can result in negative transition probabilities.